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I describe composite pulses during which the average dipolar interactions within a spin ensemble 
are controlled while realizing a global rotation. The construction method used is based on the average 
Hamiltonian theory and rely on the geometrical properties of the spin-spin dipolar interaction only. 
I present several such composite pulses robust against standard experimental defects in NRM: static 
or radio-frequency field miscalibration, fields inhomogeneities. Numerical simulations show that the 
magic sandwich pulse sequence, a pulse sequence that reverse the average dipolar field while applied, 
is plagued by defects originating from its short initial and final 7r/2 radio-frequency pulses. Using 
the magic composite pulses instead of ir/2 pulses improves the magic sandwich effect. A numerical 
test using a classical description of NMR allows to check the validity of the magic composite pulses 
and estimate their efficiency. 

I. INTRODUCTION 

Coupling of spins via the dipole-dipole interaction is of tremendous interest in solid state and liquid state nuclear 
magnetic resonance^ (NMR). It also has deep implications in the field of quantum information processing as dipole- 
dipole interactions can be used to implement quantum logic gates^ but can also produce significative coupling to a 
thermal bath leading to decoherence^. 

In solid state NMR, the important strength of the dipolar coupling between first neighbours is an important source of 
effective transverse relaxation. The "magic sandwich" (MS) invented by Rhim, Pines and Waug is a pulse sequence 
that produces an average dipolar interaction equal to minus half of the free dipolar interaction while applied allowing 
to perform so called time-reversal experiments in samples where the dipolar interactions dominate the Hamiltonian. 

The MS pulse sequence originating from solid state NMR has recently been proven useful in liquid state NMR, 
namely a hyperpolarized 3 He- 4 He liquid mixtures at low temperature 5 (~ IK). In this dense liquid, 3 He nuclear spins 
1 /2 are polarized by optical pumping up to 20% and present a strong long range dipolar interaction^ in comparison 
with solids where the dynamics is dominated by promiscuous spins. Peculiar dynamical behaviours have been obser ved, 
in particular spectral clustering at low tip angle and the instability of the uniform magnetization at large tip angle^l 
In this system, the MS has been successfully applied to stabilize the NMR evolution and should allow to gain insight 
on the development of incipient turbulent spin dynamics^. The MS induced time-reversed evolution is however 
limited in time, and numerical studies of the NMR dynamics have allowed to point to a built-in defect in the MS 
pulse sequence due to initial and final 7r/2 radio-frequency (rf) pulses which cause a significant perturbation in the 
average Hamiltonian. 

I have designed and caracterized new composite pulses-shaving the same qualitative feature as a single 7r/2 pulse 
sequence but being less intrusive on the pure dipolar evolution. Those pulses are based on geometrical constructions 
derived from the average Hamiltonian theory 11 and have the same effect as a magic sandwich pulse sequence while 
applied thus earning their name "magic composite pulses" . 

This paper is organized in two parts: First, we briefly review average Hamiltonian theory applied to the spin-spin 
dipolar interaction. We extract a geometrical sufficient condition on a pulse sequence to cause an average dipolar 
Hamiltonian proportionnal to the free dipolar interaction term. A defect of the magic sandwich due to finite pulse 
length is described and several magic composite pulses are proposed to correct it. In the second part, a numerical 
simulation algorithm allowing the computation of the NMR spin dynamics of classic magnetic moments on a square 
lattice network is decribed. Numerical simulations allow to evidence the imperfection of the standard magic sandwich. 
We show that analysis of the spin dynamics computed with this classical model allows to test qualitatively and 
quantitatively the magic composite pulses. 
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II. AVERAGE HAMILTONIAN THEORY AND CONSTRUCTION OF THE MAGIC COMPOSITE 

PULSES 

A. Average Hamiltonian 

The spin Hamiltonian in which we will be interested describes the dynamics of like-spins Ii in a static homogeneous 
magnetic field submitted to a time- varying rf field and for which we take into account the truncatecP dipole-dipole 
coupling Hd- 

Using the convention h = 1, this Hamiltonian reads 



where 



and 



H(t) = H Q + H Ti (t) + H d (1) 

H = n J2i z< h (2) 

i 

H l{ (t) = n x (t) cos (ojt + 4>{t)) 4i> ( 3 ) 

^ = EE 6 ^ Ai i ( 4 ) 

3 i<3 



with 




(5) 



and 



% = 7 ^3 ( 6 ) 

In liquid state NMR, this last quantity has an implicit time dependence due to atomic motion. 

If the spin ensemble is described by the density matrix p 1 its dynamics is governed by the Liouville equation 

^ = -i[H(t),p]. (7) 

For a static Hamiltonian, a solution of this equation can be simply computed: In this case, the density matrix at 
time tf is related to the one at time to by p{tf) = U(tf — to)p(to)U~ 1 (tf — to) where U(t) = exp (iHt) is the unitary 
propagator for the time-independent Hamiltonian H. 

For a time-dependent Hamiltonian, no simple solution of the Liouville equation exists. However, a propagator 
U(tf, to) can still be defined. It is then supposed that this propagator can be written in a form similar to the obtained 
in the case of a time-independent Hamiltonian: 

U(tf,t )=expiS(t f ,t ), (8) 

where S has the dimensionality of an actiorPM If we know the spin state at time to and are only interested in 
its state at time tf, H = S(tf,t )/(tf — t ) can be used as effective Hamiltonian instead of H(t) in the Liouville 
equation^!. 
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This expression for the propagator U may give wrong intuitions as this average Hamiltonian is only valid to compute 
the state at time tf knowing the intial state at time to- The states obtained at intermediate times have no physical 
reality and for this reason, the average Hamiltonian method is often described as a stroboscopic method. However, 
in spite of this inherent limitation it has been widely used with great success in many fields of physics. 

In NMR, the average Hamiltonian method is particularly useful as most NMR pulse sequences are periodic. In this 
case, if r is the period of the sequence, then the average Hamiltonian used for the time translation from to to to + t 
is also valid for any t n = t + nr. Using symmetry considerations, it is often possible to cancel several orders of the 
average Hamiltonian. Moreover, in the case of a periodic pulse sequence, Floquet's theorem garranties the existence 
of the form ^ for the propagator U. 

B. Magnus expansion of a time-dependent Hamiltonian 

The average Hamiltonian H is a theoretically well defined quantity but pratical determination is by no mean an 
easy task. Several methods have been described to obtain approximations of H. Among them, the most successful in 
NMR is probably the Magnus expansion^ which allows to compute successive aproximations of H ; 

H = H<® + hM+SW + --- (9) 

Various methods have been used to express the Hf"': first by Magnus^ and later by Wilcox-^ and HaerberlerP^. 
For reference the first two terms are given here considering, without lack of generality, that to = and tf = r: 



ff (o) = 1 r toxins, (io) 

T Jo 



ff (1) =^ / dta / dt2[ff(*i),ff(*2)]. (11) 



JO 



General expression for the term of order n implies n nested commutators of the time-dependent Hamiltonian H(t) 
at different times 

A sufficient condition for this series expansion to converge toward the average Hamiltonian is that the form Q 
for the propagator exists. In the case of a periodic pulse sequence existence, and thus convergence, is garrantied by 
Floquet's theorem^]. The problem for using the magnus expansion is thus rather whether truncation of the Magnus 
expansion is representative of the sum. This criterion is given by 



|JT(t)||dt<7r, (12) 



o 



where |.| is the supremum norm: 



\\A\\= sup YtT' (13) 

There are many cases however where this criterion is not met experimentally and where truncation of the Magnus 
expansion still provides accurate prediction on the experimental observations. The reason for it is that, in a large 
sample, only few modes of the average Hamiltonian are resonant with the pulse sequence spectrum, so that truncation 
fails only for these modes but most of the sample is correctly described by the truncation. In this case, the truncation 
of the Magnus expansion does not predict slow relaxation due to warming of the spins by the rf pulse sequence. This 



explains why, also in many cases assumption ( 12 ) is not verified, truncation of the Magnus expansion is used with 
success in solid state NMR. 

In the case of magic composite pulses, the pulse sequence not beeing periodic, assumption (fSp is not garrantied. 



Casas et al. have shown that criterion ( 12 ) is sufficient to have the existence of ([8j) and thus the convergence and 
validity of the truncation for the Magnus expansio n 1 17 ^ 18 !. However, in the case of a dipole-dipole Hamiltonian, criterion 



(12 1 is cxpcrimmcntally tremendously difficult to meet due to large eigenvalues in the Hamiltonian 



In fact partial validity of the Magnus expansion beyond criterion ( 12 1 is provided by the same argument used in the 
case of Magnus expansion truncation in the case of periodic Hamiltonians: In the general case, the form Q does not 
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exist for the whole system. However, it is possible to define a large subspace of the Hilbert space in which truncation 
of the Magnus expansion can be used to described the NMR evolution (only few excluded modes beeing warmed by 
the rf pulse sequence). This imply that the pulse sequence decribed in this paper have to be taken with the same 
cautiousness as regular periodic magic pulse sequences in solid state NMR as they rely on similar assumptions. 
In the case of quantum information processing, approximate validity arguments probably fail considerably and 



condition ( 12 ) must be met absolutely if magic composite pulses are used. 

For a given Hamiltonian, a change in the reference frame can change the domain of convergence of the Magnus 
expansion. A carefully chosen reference frame can also provide faster convergence of the Magnus series^. In the ca se 
of the dipolar interaction in NMR, the proper reference frame to work with is the doubly tilted rotating frameJ^^HH] 

C. The doubly tilted rotating frame 

The doubly tilted rotating frame is the frame in which the instantaneous Hamiltonian is only the dipolar 
interaction^. The (passive) rotation operator transforming the laboratory frame in the doubly tilted rotating frame 
is the propagator computed by taking into account the external and applied magnetic field only. It can be deduced 
from the Liouville equation and satisfies the equation : 

id t R = (H + H Ti {t))R (14) 

with i?(0) = Id (the instant t = is chosen appropriately). The Hamiltonian in (JtJ) can be replaced by H'(t) in the 
doubly tilted rotating frame where H' '(t) is given by 

H' = R(t)HR- l (t) - idtR^R- 1 ^). (15) 

Using this transformation on the spin Hamiltonian ([!]), we get 

H'(t) = H' d {t) = R(t)H d R-\t), (16) 

as extra terms corresponding to the applied field have been canceled by going into the doubly tilted rotating frame. 
We replace Hd by its expression given in equation Q and apply the following simplification : 

RtliALjR- 1 = {Rt^R-^iRAR-^iRfjR- 1 ) (17) 

at this point we remark that R(t) is an evolution operator which is composed of infinitesimal rotations in the SU(2) 
space and is consequently a rotation operator this space. We call R the corresponding rotation in SO(3). So that we 
get 

RljR- 1 = (Rfj), (18) 

and 

Rtf-R- 1 = *(r8). (19) 
Moreover since R alone is a scalar quantity in the SU(2) space : 

RAR' 1 = A, (20) 

hence 



with 



(21) 
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1 

A'(t) = R 1 ( 1 | R. 

0-2 



(22) 



Finally, we obtain the following expression for H'At) which we will use whenever we apply the Magnus expansion: 



^(<) = EE^A'(«)/] 



(23) 



J i<3 



Let us remark that equation (22 1 implies that the matrix A!(t) can be written in the following simple form 



A'(t) = Id-3*n(i) -n(t) 



(24) 



where n(t) is the unitary vector pointing to the z direction of the laboratory frame from the doubly tilted rotating 
frame. It is thus solely determine by 2 angles. 



D. Magic conditions 



The computation of H 1 ^ is a priori complicated, but we remark that the sequence operator R(t) only concerns 
the vectorial part of the Hamiltonian (U| and does nothing to its spatial dependencies The term 3'} ^ s then 
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1(0) 



EE* 

j i<3 



(25) 



We define a magic sequence as a sequence which rotates the spins and does not modify the dipolar interaction 
during a given duration (except by a scaling factor). The zero order condition for magic sequences formally reads 



H 



'(o) _ 



k3. 



(26) 



It can easily be shown that — 1/2 < k < 1. 

From this condition we can extract a general geometrical condition by separating the spatial part from the geomet- 
rical part of the dipolar interaction. We get the following geometrical condition: 



- [ A'(t 1 )di 1 = fcA 

T Jo 



(27) 



It can be simply verified that conditions ( 26 1 and ( 27 ) are equivalent 



For higher orders the condition of order n for a magic sequence is simply: 

3'} n) =0 if n > 1 



(28) 



As for order 0, we show in appendix [X] that geometrical conditions on the pulse sequence exist and are sufficient 
to ensure that (28) is verified. The equivalence of those geometrical conditions with the corresponding conditions on 



the average dipolar spin Hamiltonian is not prooven in the general case (although appendix |E] equivalence for orders 
and 1 is prooven). 



E. Taking into account experimental defects 



Experimental NMR pulse sequences are plagued with defects. The idea to deal with it consists in separating the 
time-dependent Hamiltonian ([T]) in a ideal part and a perturbative part : 
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H{t) = H + flrf(t) + H d + SH + SH lf {t) (29) 

Here H is the ideal static field being equal to toI z and 5H is the difference with the applied field. Similarly H T f(t) 
designates the ideal pulse sequence following the theoretician expectations and 5H r f(t), the errors in the experimental 



realization. We now express Hamiltonian ( 29 ) in the doubly tilted rotating frame 



H'(t) = H' d {t) + R^SHoR' 1 ^) + RtySHrftyBr 1 ® (30) 

Applying the average Hamiltonian theory we get an intricate average Hamiltonian H' . However at the lowest order, 
the average Hamiltonian is a simple sum of the three different effects: 

H>W=H'W+5H'W+6Hf ) (31) 

1. Robustness against a Larmor shift 

If the local Larmor frequency is not equal to the rf, the perturbation brought on the spin dynamics has to be taken 
i nto account. The extra term in the average Hamiltonian due to a Larmor shift can be easily calculated using equation 

|lo| with the Zeeman term 5Hq = Sflo I ' z. We obtain the effective hamiltonian term : 

5H ' = -[ R(t)SH R(t)- 1 dt=- f dt 5tt (R(ti)(f ■ z)R(t 1 y 1 ) (32) 
T Jo T Jo 

We use the same simplification technique : 

R{T- z)R- 1 = {RIR- 1 ) ■ z = (RJ) ■ z = T ■ (R^z) (33) 

This last expression allows to write a simple geometrical criterion for a sequence to be robust against a Larmor 
shift: 



R" 1 (t)zdt = (34) 



This criterion can be verified easily for a given pulse sequence. In NMR experiments, composite pulses robust 
against a Larmor shift are desirable to prevent failure of the pulse sequence due to experimental imperfections of the 
field homogeneity or the application of magnetic field gradients. 

2. Robustness against perturbation in the excitation field B\ 
In the case of a perturbation in the excitation field Bi, the condition of robustness at the lowest order is 



1 r T 

5~H[ ( { 0) = ~ R(t) &Hrf (*) -R(f ) ~ 1 df = 0. (35) 
T Jo 

Unfortunately in this case no obvious geometrical equivalent condition exists because of the time dependence of the 
perturbative term 6H r f(t) which depends on the pulse sequence used. Moreover, if solutions exist for a miscalibration 
of B\ amplitude, no possibility can be found in the case of a miscalibration of phases (except for pulse sequences 
doing a 2tt rotation). 
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F. Examples of magic sequences of order zero 

1. The magic sandwich 

The MS^J consists in applying a strong rf between to perfect tt/2 pulses. Here we compute the effective Hamiltonian 
ff(o) during the continuous rf field along the x-axis. The rotation operator R x (i) 

/ ° ° \ 

R x (f) = cos^ sin^ , (36) 



-sin%^ cos^ 



so 



A'(t) = I -i + fcosQxi fsinftii ). (37) 



3 sinfM -i-fcosfii* 
If we do a n rotation the average interaction matrix A' is 




A' = -I . (38) 



A' is not directly proportional to A, but is with k = —1/2 when sandwiched between two tt/2 rotations along the 
y-axis (this structure is the reason for naming this pulse sequence "magic sandwich"). 

2. Effect of the n/2 pulses finite duration in the MS 

If we suppose that pulses are not infinitely short, we have to take into account the imperfections brought by the 
initial and final tt/2. The initial tt/2 consists of a rotation along the y-axis. We compute the effective interaction (the 
pulse last t = 7r/f2i) : 



A'{t)dt= 10 (39) 




o 

This interaction matrix is not proportional to the free interaction matrix because of those perturbative non-diagonal 
terms. Consequently, the first and last n/2 pulses finite duration of the magic sandwich are sources of perturbation 
on the ideal NMR evolution. An illustration of the effect of those pulses by numerical simulation is given in section 

Unci 

3. A 7r/2 magic composite pulse 

In order to suppress the perturbation introduced by the tt/2 pulses in the MS sequence, we replace the simple 7r/2 
rotation by the sequence 

(tt* - (tt/2), - (7r/2) V) i/a - (7r/2)_ y , 1/2 ). (40) 

The indexes 1/2 indicate that the pulses durations are the half duration of the standard pulses. The total duration 
of rf pulses along x is equal to the total duration of rf pulses along y. We compute the average interaction matrix by 
dividing the computation in two terms: the first part consisting of a a; rotation and the second part consisting of a y 
rotation : 

A' = 1(A x + a;) (4i) 



The geometrical term corresponding to the tt x pulse is given by equation ( 38 1 . The rf along the y axis is more 
complicated, but it has been shaped so that the time spent in every direction on the circle described on the Bloch 
sphere is the same and consequently at order 0, causes the effect of a 2ir homogeneous rotation^. 

-\ 

A; = I 1 1 (42) 



-\ 



This concludes the calculation on the effect of the magic pulse sequence: 




A' = | i | (43) 

the effective dipolar interaction at the order zero is k = 1/4 and the equivalent rotation during the pulse is n x 
- (it /2) y which brings the z component of the spin along the x direction. This property makes this pulse sequence 
suitable to replace tt/2 initial and final pulses in a magic sandwich. 

4- Other magic composite pulses n/2 

Using the technique described in the previous section, we can build other Mir/2 which can be more efficient in 
certain situations. I give few of them here (proofs are left to the reader) 
The first inhomogeneous rotation along the x-axis can be chosen between: 

2tt x , 2tt 

— X 

2tt x ,2tt- x (44) 
e 2^±x 

The second inhomogeneous rotation along the y-axis can be chosen between the following sequences: 

7r/2-y,i/2,7r/2y,i/2,7r/2 y 

ir/2 y ,TT/2 yA/2 ,TT/2_ yA/2 (45) 
a-j/,1/2) a y,i/2, K/2 y , (tt/2 - a) V) i/ 2 , (tt/2 - a)^ y l/2 

Those last three pulse sequences have been shaped to give an average dipolar matrix equal to the one obtained 
using a single 2~k v pulse. By assembling two inhomogeneous rotations among (44) and (45) having the same duration 
along the i-axis and y-axis, a magic composite tt/2 pulse is formed. 

Several variations are possible, since we only chose one particular strategy. For example, we could decide to mix x 
and y rotations in a complex way, or use a free evolution period to compensate for diagonal terms. Here is an example 
where the second part of the sequence is only composed of pulses of the same intensity: 

Kx,7/6 > n/6-y , tt/6 v ,Tr/2 y , tt/6 v ,ir/6^ y (46) 

This sequence has the advantage of being shorter (for the same duration reference) than the first Mir/ 2 and is then 
expected to be more efficient. This time, the y sequence is not a simple equivalent to a 2n homogeneous rotation, 
although it has the same properties (see appendix [B]) . 



Other magic composite pulses robust against a Larmor shift can be created by using the extra condition ( 34 ) . The 
following sequence is both magic and robust against a Larmor shift and realize a total rotation equivalent to tt /2- y 
pulse: 

2tt x , 3tt/2 v , 7r/2.y :0 .5, ~^/^y,o.5 (47) 

All the magic composite tt/2 pulses presented until now are of order 0. In appendix [C] I describe the construction 
of a Mir/ 2 of order 1 : 
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[(lT, ~7r)y :2 .767, (l", -1")x,2.767, ("", -7r)x,2.767] j [%,1.488j (l"/4l.557 > 7r/4l. , 7r/4 2 .557) a , 

(tt/42.557, 7r/4i. , 7r/4i.5 57 ) a , 7r Vj i.488] 

This pulse sequence is particularly long and difficult to realize with an NMR spectrometer but will nontheless be 
used in numerical simulations to illustrate the efficiency of a numerical test for magic composite pulses. 



G. Enhanced magic sandwich 



Using magic composite 7r/2 it is possible to enhance the efficiency of the MS (see section III C ). The usual MS pulse 
sequence is 



7r/2 3/ ,(2 7 r :c ,-27r x ) n ,7r/2 3/ , (49) 

where the filling with (27r, —2tt x ) pulse sequence has been chosen to avoid long term dephasing effects due to a 
miscalibration of B\ amplitude and a local error in the Larmor frequenc}^. 

To obtain an efficient sequence, the normal 7r/2 rotations are replaced by magic pulses. We call the resulting 
sequence an Enhanced Magic Sandwich (EMS), since it should be more efficient in reversing the average dipolar 
interaction without introducing perturbations: 

[the, tt/2^, 71-/2^0.5, — 7r/2 yj0 . 5 ], (277^, — 2-K x ) n , [n/2 y , 7r/2 yi0 .5, — tt/2 

y,0.5 •> f^x 

} (50) 

1. Larmor robust zero order Mix 

M7T pulse sequences are interesting as they allow population inversion. The construction of such pulse is easier 
because less symmetries are broken by the total rotation needed which allows more flexibility in the design. In 
particular using only 7r pulses of phases <pi, 02, etc., the geometrical criterion for magic condition of order is simply 

^>e 2 ^=0, (51) 

t 

where Ti is the duration of the zth w pulse. Likewise, the criterion for Larmor shift robustness of order is 



5>e^=0. (52) 

i 

We deduce several M7r of interest for usage in experiments, for example in multiecho CP-like pulse sequence d 23 * 24 ! 
where the dipolar interaction should not be sensitive to the pulse sequence (appart from a scaling factor). 



Pulse sequence 


Magic 


B 


B 1 




X 


X 


X 


277^, TT y , TTyfi.5, — %, 0.5 


/ 


/ 


X 


n iT/6l n TT/2y 7T 5TT/6 


/ 


/ 


X 


K-Tr/61 1t tt/27 7t -5-!t/6 


/ 


X 


/ 


2tT-tt/2, 7r 7r/2, ^5tv/6 


X 


/ 


/ 


2tT-tt/2, TT^/2, 77^/6, 77^, 7r 57r / 6 , 7T 


/ 


/ 


/ 



Magi c com posite 7r pulses have been reported in the litterature as allowing population inversion in dipolar coupled 
systemd 25 * 26 ! . In particular, the pulse sequence Tr.^/g, 77^/2, 77_ 57r / 6 have been studied by Tycko both theortically and 
experimentally^!!. 
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III. NUMERICAL TEST OF THE MAGIC COMPOSITE PULSES 

In this section, I first explain why a numerical simulation of dipolar coupled magnetic moments is relevant to test 
the efficiency of magic composite pulses ususally applied in the context of quantum mechanics formalism. I will then 
briefly describe the simulation algorithm and present numerical evidence for the perturbation brought by the MS 
on the theoretical NMR evolution. Finally, I will describe the design of a numerical test for magic composite pulses 
and present the results on several of them. This test is a useful tool to compare the efficiency of the different magic 
composite pulses. 

A. Equivalence between the quantum and classical description of NMR for the design of magic sequences 

The nature of the dipole-dipole interaction naturally leads to the formation of entangled states and for this reason, 
a classical description of NMR seems inappropriate. In liquids however, an hydrodynamic approach is suitable to 
describe the NMR observations because of the fast relative motion of promiscious spins and the consequent loss of 
quantum information^. Both experimental confirmations of predictions from the classic formalism and theoretical 
advances now support that the classical formalism correctly describe the NMR evolution. In this formalism, the spin 
ensemble is described by the local magnetization m(r, t) whose evolution can be described in the rotating frame by 
the mesoscopic Bloch-Redfield equation^ which takes into account the coupling between magnetic moments due to 
the distant dipolar field: 

d t rfi{r, t) = 7 m(f, t) x B d (r, t) (53) 

with 

Mr,t) = £ J^^-^kl l)Am(rV) (54) 

where k = pErpq" an d the volume of integration V is the total space M 3 minus an infinitesimal spherical volume 

centered on the point r. The distant dipolar field is meaningful in liquids due to motional averaging of the short 
range dipolar interaction. This dipolar field brings a non linear and non local local term to the Bloch equation. In 



general, no simple solution of (53) exists for the dynamics of the magnetization. 

The Bloch-Redfield equation gives the evolution of the classical description of NMR involving coupled magnetic 
moments. In essence, it is truly different from the quantum description of NMR described earlier. In liquids however, 
this formalism describes accurately the NMR evolution. 

We arc interested in using the magic composite pulses described in the previous part. Despite the different nature of 
the quantum ([7| and classical (53) NMR equations of evolution, we show in appendix [d] that the sufficient geometrical 
conditions presented earlier in the quantum formalism are still valid for classical magic conditions. Furthermore, for 
orders and 1, quantum and classical magic conditions arc exactly equivalent as shown in appcndixjE] This equivalence 
between the classical and quantum magic conditions ensure that a numerical simulation of dipolar coupled magnetic 
moments can give relevant information about magic composite pulses used in a quantum formalism context for example 
for quantum information processing. Using this property we present below a numeric test of the magic composite 
pulses. 

B. Numerical Simulation 

To compute NMR dynamics, a general-purpose model is used based on an exact calculation of the time evolution of 
coupled magnetic moments on a cubic lattice. This enables to account for all the terms acting on in in experimental 
conditions in an exact manner, including diffusion, dipolar interactions, rf field, and coupling to the detection coil. 
However, this coarse-grained description of a continuous fluid has obvious limitations when spatial variations at the 
scale of the lattice constant are involved. Computations are performed in the Larmor rotating frame using a standard 
secular approximation. The magnetic field induced at each site by the remainder of the sample is efficiently computed 
by toggling between real and Fourier spaces using a half complex FFTW transforrrPS. The time evolution is computed 
by integrating the Bloch-Redfield equation using a Runge-Kutta 4-6 algorithm modified to be robust against evolution 
discontinuities (due to rf pulses). Free evolution usually involves slow local changes and computations are quite fast 
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FIG. 1: Numerical simulations of the NMR evolution of spins submitted to MS pulse sequences. The signe is not given during 
the application of the MS as experimentally the rf pulse sequence saturates the detection circuit. An ideal time-reversal is 
given to allow a comparison of the pulse duration influence on the observed echo shape. An equivalent amplitude of 1000 turn 
per second during the MS filling is sufficient to reduce the perturbation brought by the filling to a negligible part of the total 
perturbation. The delay and amplitude of the echo is directly related to the finite duration of the initial and final n/2 rotations 
indicating that these are the cause for the perturbation of the evolution with regards to the ideal time-reversed evolution. 
Consequently, an EMS (see section II G I improves the echo shape while using a lower equivalent amplitude during the initial 
and final 7r/2 rotations. 



on standard PC computer (Core2Duo 6700). When intense rf fields are applied, the computing load is significantly 
increased (e.g. by one order of magnitude for MS sequences). The computations described below have all been done 
on lattices of 32 x 32 x 32 sites. Because of periodic boudary conditions, only 20 x 20 x 20 sites were containing 
magnetic moments in order to account for finite size effects. 



C. Numerical evidence for a defect in the magic sandwich 



In the presence of long-range dipolar spin interactions, the spins in a liquid can present various collective behaviour. 
After a 7r/2 tipping pulse, the distant dipolar field (DDF) produces NMR precession instabilities which cause the signal 
to collapse. A model developed by Jeener and numerical simulations have shed light on magnetization dynamics^ 
at large DDF, during transverse precession, unstable inhomogeneous magnetization patterns develop and grow expo- 
nentially in time at a rate T. These patterns arise from an initial seed of inhomogeneity following the tipping pulse. 
This instability produces a brutal decay of the signal. It has been showrP that applying a magic sandwich in this 
liquid system allows to produce time-reversal evolution and generates a dipolar echo. However this echo of the initial 
signal is imperfect being both damped and delayed. Figure [T] represents numerical simulations of this experimental 
situation. Different MS have been used including an EMS using magic composite pulses for initial and final 7r/2 
rotations. Those simulations reveal that the echo imperfection is mainly the result of the initial and final 7r/2 finite 
duration (continuous rf applied during the MS having little effect in experimental conditions). It also shows that 
using a magic composite pulse of order zero is equivalent in the perfect conditions of the simulation to reduce the 7r/2 
pulse duration by one order of magnitude. 



D. Design of a numerical test for magic composite pulses 

I consider the following numerical experiment (see chronogram in Fig. |2j: A free evolution followed by the pulse 
to be tested. Right after the pulse, an instantaneous rotation compensates its global rotation and is followed by a 
time- reversed evolution at half the pace of the free evolution (in order to mimic the effect of a magic sandwich). 
The perturbation introduced by the use of the tested pulse sequence instead of a perfect instantaneous rotation is 
evolving in an unknown way as long as the magnetization is mainly inhomogeneous. Around the instant of the 
expected echo, the signal amplitude is closed to its initial value, implying that magnetization is homogeneous agairP^l. 
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Time (s) 

FIG. 2: Chronogram of the numerical test for magic composite pulse and example on a M 

pi/2. The different black curves are obtained for different magic composite pulse durations (durations are divided by two 
between two consecutive curves ranging from 0.5 s to 0.015625 ms). The deviation observed between the perfect evolution 
(dashed line) and the numerical tests for finite pulse sequences is due to the presence of a parasitic inhomogeneity introduced 
by the composite pulse. The fall time measured at a reference amplitude of 90% is a growing affine function in the logarithm 
of the Bi intensity. The dashed square indicates the valid zone where Tf measurements are relevant and is defined by the 
conditions given in the text. 



For a sufficiently small initial perturbation brought by the tested pulse, the perturbation remains small during the 
evolution pre ceeding the echo instant and, as a small perturbation, evolves according to a perturbative treatement of 



equation ( 53 ) in a linear manner driven by the main magnetization map^l. When those two conditions are reached (in 
numerical an initial free evolution time of 300 ms for 30 Hz of characteristic DDF ensures that they are), the remaining 
magnetization perturbation grows exponentially during the echo as it does during the intial phase of evolution. By 
measuring the collapse of the signal during the echo, we thus obtain quantitative information about the perturbation 
of the magnetization map brought by the tested pulse sequence. 



In figure [TJ the numerical is applied to the M7r/2 described in II F 3 for various pulse durations. We note the fall 



time Tf when the signal cross a given amplitude after the top of the echo. The relevance of this measurement is 
determined by several conditions: 

• The echo amplitude should be at least 99% of the initial magnetization implying a mainly uniform magnetization 
during the echo; 

• The reference amplitude for measurement shall not be too low to ensure the perturbation brought by the pulse 
evolves in a linear manner. During the tests, 90% appeared to be an acceptable value; 

• The echo does not reproduce the ideal echo shape as, in this case, the fall time is not related anymore to the 
perturbation introduced by the pulse. 

For this window of parameters, the fall time Tf represents a global measurement of the efficiency of the pulse sequence: 
A benchmark of known magic pulse sequences allows to obtain numerically the order of the sequence in experimental 
condition. Examples given below show that this measurement of the sequence magic order is both precise and robust. 
This test also allows to compare the efficiency of pulse sequences having the same magic order. (This numerical test 
simulating the NMR dynamics during a magic sandiweh.) 



1. Extracting informations from the measured fall time 



During the echo, magnetization beeing nearly homogeneous, the exponential development of the magnetization 
instability is the dominant dynamical process explaining the observed echo shape. Hence the fall time Tf for a floor 
of 90% is given by the following equation: 



1 - aexpTTf = 0.9 



(55) 
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Where a is proportional to the size of the perturbation introduced by the tested pulse sequence. We thus have the 
relationship 

logO.l-lna 

Tf = p (56) 

which imply that Tf is a linear function of log a. For a magic pulse sequence of order n and duration t, the 
inhomogeneity size a is proportional to r n+2 , the first non canceled perturbative term in the Picard iteration of the 
Shroedinger equation involving n + 2 dipolar Hamiltonians. Expliciting this dependency in the sequence order we get 
the relationship 

T f oc lna(r ) - (n + 2) In t/tq. (57) 

By doing several numerical simulations for different pulse duration, it is possible to measure the magic order of the 
tested pulse. The duration r is not necessarely the duration of the sequence but may be set as the duration of a 
typical 2-7T reference pulse. By fitting Tf as an affine function of ln(r): 

Tf — A + B ln(r), (58) 

we get a measurement of the sequence's order trough coefficient B and of the relative sequence efficiency trough 
coefficient A. 



2. Completeness of the test: spin turbulence 



This numerical test using a classical formalism may seem incomplete to test the ensemble of geometrical conditions 
required for a magic pulse sequence of order n. In appendix [Ej we show that the magic geometrical conditions are 
equivalent to the magic classical conditions, the question of completeness of the test then lies in the ability of the 
numerical simulation to test every possibility offered by the magic geometrical conditions using a certain realization 
of magnetization at the instant when the pulse is applied. Although it is not possible to directly proove this point as 
the magnetization map is not controlled during the free evolution, it is highly unprobable that a geometrical condition 
would not be tested by this numerical simulation. This results from the onset of spin turbulence that magnetization 
undergoes before the test which corresponds to the formation of a quasi random magnetization map. Hence, this 
map imply many random classical conditions, and as stressed out in appendix [E] linear combinations of those classic 
conditions allow to recover the geometrical conditions. Given the huge number of random classical conditions tested, 
the probability of missing a magic geometrical condition is thus of measure zero. Of course symmetries of the system 
could default this argument, but numerical simulations are realized with an initial inhomogeneity seed chosen to break 
every possible symmetry of the system to avoid such problem. Consequently, the classical numerical test completely 
tests the magic geometrical condition and can be used safely. 

E. Benchmark of magic composite pulses 

The three following pulse sequences will be used as references for the magic order: 



Pulse sequence 


Magic order 


n/2 x 


-1 







T^xi 7Tyi ^T— yi TT— x 


1 



The magic order —1 given for the pulse sequence ir/2 x exists because of the particular choice of order labelling in 
the average Hamiltonian theory. This order has to be interpreted as giving a perturbation scaling as r n+2 = r. Figure 
[3] represents the fall time Tf as a function of the pulse equivalent duration. 

On Fig. [4jthe slopes B are represented as a function of the theoretical order n and will serve as references to define 
the measured order of the magic composite pulses. A fourth sequence has been added on this figure: a tt rotation. 
This pulse presents more symmetries than the tt/2 rotation an is thus of magic order in duration but produces less 
inhomogeneity. For this reason, the measured slope of the sequence is more important than expected theoretically. 
This pulse serves as a guide to estimate what are the maximum possible discrepancies to be expected between theory 
and simulations. 

Let us now turn to the following tested magic composite pulses: 
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FIG. 3: Fall time as a function of the duration of a 2n equivalent pulse for the three reference pulse sequences n/2 x , iT x ,iTy, 
and iT x , n y ,-K- x ,n- y . The numerical test is realized for 300 ms of free evolution and a caracteristic dipolar frequency of 30 Hz. 
A small but relevant correction was applied prior to representation to take into account the effective time-shift on the echo 
instant due to the finite duration of the tested sequence. (A correction which is 3 times the duration of the tested composite 
pulse.) 



200-, 




-1 1 

Pulse sequence theoretical order 



FIG. 4: B coefficients of the reference composite pulses measured from the data of Fig. [3] function of their theoretical magic 
order. As expected, the slopes scale linearly with n + 2 where n is the theoretical magic order of the magic composite pulse (see 
eq. (57l). A fourth pulse n x is also studied to infer the maximum discrepancies between the test results and the theoretical 
expectations. This benchmark allows to test numerically the magic order of a given composite pulse. 



Name 




tt/2 




Mfa 


Kx, - 


Mfb 




Mfc 




Mfd 


7: x,l/2, — n x 


Mfe 


7^,7/6! _7r / 


Mf 1 





Pulse sequence 



tt/2 



-tt/ 2 j/,i/2> n/2 y ,i/2, it/2 y 
■/2 y ,Tr/2 y l / 2 , -Tr/2 yl/2 



see section 



C2 



Order 



-1 






1 



Figure [5] represents the measured parameters A and B for those pulse sequences. 



F. Discussion of the efficiency of the different magic composite pulses 



The scale of the perturbation size brought by a magic composite pulse scales with its duration, in experimental 
situation, a shorter pulse will in general give a smaller perturbation. A short pulse sequence is however not always 



15 



880 
860 
840 
820 
800 
780 
760 
740 
720 
700 



Order -1 



Order 



Order 1 



M|a-d 



Mfe 



M| order 1 ■ 



20 40 60 80 100 120 140 160 
B coefficient (ms/l(ms)) 



FIG. 5: Measurement of coefficients A and B of the different tested magic composite pulses. The numerical test is realized 
in the same conidtions as for the reference benchmark. Every simulation is realized by defining a relevant equivalent scale 
of duration for the 27r pulses. Vertical lines represents the expected slopes B for the different theoretical orders —1, and 1 
extracted from Fig. [4] For a given magic order, the coefficient A quantify the relative efficiency of the magic composite pulses 
for a given 27r reference duration. 




"oi ' 1 
27t pulse duration (ms) 



FIG. 6: Comparison of the fall times function of the equivalent 27r duration for the M7r/2a and the M7r/2 of order 1. At low 
rf intensity, i.e., long 27r equivalent duration, the zero order sequence is more efficient than the first order sequence. 



desirable for several reasons. For example in low temperature experiments, the amount of rf power can warm up the 
studied sample. Moreover, in the limit of ultra-short pulses (typically Qi/Qq > 1/100), rotating wave approximation 
is not correct anymore and corrections on the NMR dynamics have to be taken into account. Due to those corrections, 
the effective order of the pulse sequence may be brought back to — 1; for example, for ultra short pulses a plain tt 
pulse does not produce an homogeneous rotation anymore. 

We have seen that the numerical test allows to determine reliably the pulse sequence magic order, it has been used 
to verify that the construction of magic composite pulses was correct. The similarity of the test with the evolution 
expected during a MS is of equal importance as the coefficient A related to the size of the perturbation for a reference 
2tt pulse duration allows the comparison between different pulses in a situation close to the experimental situation: 
The echo observed in the numerical test corresponds to the echo that forms during the MS. Comparisons are relevant 
between pulses of the same order where the efficiency is determined by coefficent A (the smaller the better) , but it is 
also relevant when choosing between pulse sequences of different orders. Indeed, depending on the duration the pulse 
sequence having the highest magic order is not necesseraly the most efficient for a given equivalent 2n duration. Figure 
[2] illustrates this and show that the magic composite pulse of order 1 described in appendix [C] is not more interesting 
than a magic pulse sequence of order in the experimental situation^ where the 2ir duration is fixed between 2 and 
4 ms. This numerical test is a method of choice to determine which, between several magic pulses candidates, is 
susceptible to be the most efficient magic pulse in given experimental conditions. 
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IV. CONCLUSION 



In this paper, we have diagnosed that the magic sandwich is limited by the finite duration of its 7r/2 initial 
and final pulses which introduce perturbations in the time-reversal evolution. We have shown that, using the average 
Hamiltonian theory, it is possible to find sufficient geometrical conditions for composite n/2 which, if respected, greatly 
reduce the perturbation introduced when the pulse sequences are applied. We have presented a construction method 
for those magic composite pulses which, as opposed to the usual strategy in NMR, does not rely on the symmetries to 
cancel those conditions but on a direct calculation of the geometrical conditions. We then have presented a numerical 
test involving the numerical simulation of dipolar coupled magnetic moments on a lattice. After explaining how such 
test, though involving a classical formalism, is relevant for the study of magic composite pulses in a general context, 
we have explained how to obtain quantitative information about the pulse sequence efficiency. The results of the 
work presented here will be useful for solid state NMR and for the study in liquid state NMR of the complex effect of 
the distant dipolar field in high resolution spectrometers. Beyond NMR, those results may be of particular interest 
in the context of quantum information processing where dipolar couplings could be controlled using coherent control 
techniques similar to NMR to allow both deterministic coupling between qubits to implement quantum logic gates 
and decoupling of the relaxation channel between qubits and a thermal bath coupled by dipolar interactions. 
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Appendix A: Existence of geometrical sufficient conditions for magic pulses of order n 



The average Hamiltonian of order n, H ( n ) , obtained by Magnus expansion^ imply a sum of n + 1 time integrals 



of n + 1 time-dependent Hamiltonians in n nested commutators. Using expression (23) for the truncated dipolar 
Hamiltonian, it is possible to split this expression in two parts: 



• A spatial part involving the coefficients bij and the spin operators I a ; 



A geometrical part involving the matrices A(ti) and the n + 1 time integrals. 



As the geometrical part depending solely on the pulse sequence and not on the sample studied, it implies that cancelling 
this geometrical tensor Tki,k 2 .- - .k 2 „+2 (^i = x,y or z) is a sufficient condition to ensure that H^ n ' — 0. Application 
to order 1 is presented in appendix [C] It has to be noted that many conditions on the tensor T are redundant due to 
symmetries of the dipolar interaction and the average Hamiltonian expression (see for example appendix [E| . 



Appendix B: A useful geometrical transformation for magic composite pulse construction 



The computation of the quantity 



A' = -/ dt R" 1 (6»(i))AR(6»(i)) (Bl) 



T 



can be simplified by defining 



V (6) = -J- (B2) 

^ ' T Q{t) ^ ' 



and make a change in the integral variable: 



A'= / d6 r)(6) K- 1 (9)AR(9) (B3) 
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where 8m and 9 m are the maximum and minimum angles reached in the Bloch sphere during the time interval. 

For a usual rotation such as a tt x for which the rotation is homogeneous and hence rj(9) = 770, this average interaction 
matrix is easy to compute. At the magic order 0, the weak (see appendix[C]) magic condition is fulfilled if the barycenter 
of the angular distribution density rj{9) is null. This representation gives a helpful picture of the sequence effect. 

Appendix C: Construction of a magic pulse of order one 
1. First order magic geometrical condition 



Recalling equations (111 and (16), the magic condition of order 1 is 

H'P = ? f d *i P [H' d (h),H' d (t 2 )] = 0. (CI) 



2r Jo 



n 



This condition is automatically fulfilled if the pulse sequence is time-symmetric^, but magic composite pulses 
realizing a global rotation cannot possess such symmetry. Using the technique described in appendix [XJ we derive 
the magic geometrical conditions of order 1: 



For every (a. b, c, d) 



/ dti f dt 2 (A' aib (t 1 )A' c , d (t 2 ) - A , a , b (t 2 )A' C)d (i 1 )) = 
Jo Jo 



2. Construction of a first order magic tv/2 



(C2) 



To construct the pulse sequence, we will use a recursive strategy: It consists in using the magic pulse of order zero 
and modify it until it verifies the first order too. This way, the zero order conditions are already fulfilled (which is a 
great simplification in the magic conditions of higher order) . I make use of the symmetries in the construction process 
to increase the efficiency of the search by automatically fulfilling some conditions. 

As for magic composite pulses of order zero, we decompose the construction process in inhomogeneous rotations 
along axis x and y. For an inhomogeneous rotation, the magic conditions of order zero are not well defined, but we 
start with a weaker condition beeing that the average matrix A' is diagonal. If we call 9(t) the angle of the rotation 
function of time, we have the decomposition for A'(t): 

A'(i) = A' + cos0(t)C + sin(9(i)S (C3) 

and A' is diagonal if and only if 

/ cos9(t)dt = (C4) 
Jo 

[ sin0(i)di = O. (C5) 
Jo 

We will call those conditions weak magic conditions of order 0. Supposing those two conditions are fulfilled and 



introducing (C3) in (C2), we get three geometrical conditions on the rotation to fulfill the first order magic condition 
for the inhomogeneous rotation: 

/ tcos9{t)dt = 0, (C6) 
Jo 



[ tsm9(t)dt = 0, (C7) 
Jo 
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FIG. 7: Schematic of the first order magic pulse construction process. Conditions in black are fulfilled by solving the 
corresponding equations while grey [online red] conditions are solved by using the symmetries of the pulse sequence. Diagonal 
conditions of order and 1 are the remaining conditions implying necessarily rotations in both x and y axis allowing to turn a 
weak magic condition in a full magic condition. 
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and 



dii / dt 2 sin(e(t 1 )-fl(t a )) = 0. 



(C8) 

Despite simple expressions, it is not so easy to find solutions to this set of equations. We define that a sequence 
is of the first order in a weak sense if the tensor coefficients involving at least a non diagonal term of the matrix A' 
are all null. We then use the following building process: We start by building two kinds of inhomogeneous rotations 
that are O-weak/l-weak: We use composite pulses realizing 7r/2 rotations and 2tt rotations (a simple %, — 7r pulse pair 
respect the given criterion for the latter case). We then bring together a 2-k inhomogeneous rotation along x and a 7r/2 
inhomogeneous rotation along y having both the same total duration to obtain a magic composite pulse of order 0. 
As mentionned earlier, time-symmetric pulse-sequences automatically verify first order magic conditions. As we need 
a global rotation, we cannot directly use symmetrization here, but this symmetry argument is valid in a weak version 
if first order weak magic conditions are fulfilled (which is the case here). Hence we use two extra inhomogeneous 
rotations before the previously built pulse sequence to obtain a complete magic composite pulse of the first order. 

The building process is described graphically on Fig. [7j The difficult step consists in building a O-weak/l-weak 
inhomogeneous rotation. I have realized this by using 7r/4 pulses of different durations and using a time-symmetry to 
automatically fulfill two magic conditions. In practice, once the the pulse structure is determined, non linear equations 
on the timings are obtained from the magic conditions (C4), (C5), (C6), (C7), (C8) that take the form of zeroes of 
polynomials. It is clear, that without symmetry arguments, at least 5 different durations must enter in the 5 magic 
conditions in order to fulfill all of them. By using symmetries of the pulse sequence, I have been able to reduce this 
number to 3 independent parameters keeping the polynomial equations to the second degree. 

I obtain the following magic composite pulse: 



[(•7T, -7r) y ,2.767, ("", —?t)x,2.767, -71") 



a;,2.767J 



[%,1.488 5 (7r/4l.557,7r/4l., 7T/4 2 .557) a , 
( 7r / 4 2.557, T 1 "/ 4 !. , 7r/4i.5 57 )„ , TTy^ASs] 



(C9) 



Where the timing coefficients are zeroes of polynomials determined by conditions (|C4j), (C5), (C6), (C7), (C8|. 

The first order magic 7r/2 described above is a complicated sequence. Moreover it is experimentally useless, but 
it allows a convincing demonstration of the geometrical construction method reported in this paper. This construc- 
tion process also illustrates how the necessary drop of the symmetric construction strategy yields to more complex 
construction methods. 



Appendix D: Picard iteration and magic conditions in the classic case 



Working with the Bloch-Redfield equation rises the problem of dealing with a non linear differential equation instead 
of a differential linear equation in the quantum case, hence Magnus expansion is not usable. The specificity of magic 
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conditions i. e. the cancelling of the average Hamiltonians for order greater than allows to use the identification of 
the Picard iteration of the differential equation instead of Magnus expansion to obtain the magic conditions. Starting 
from a non linear differential equation of the kind: 



d t Y = A{Y,t)Y, (Dl) 

we can suppose that an operator A T (Y) exists such that the solution of the non linear equation d t Y = A T (Y)Y 
gives a solution Y(t) equal to the one obtained with (Dl|. Identification of the terms in the Picard iterations are 
possible if A(Y, t) is a linear function of Y. In the case of the dipolar interaction, this requirement is verified. If we 

suppose that magic conditions up to order n are fulfilled, this mean that A T (Y) = Ar(Y) + 0({Ar) n+1 ) where A 
quantify the size of the matrix A(Y, t). Introducing this hypothesis in the Picard iteration, we get the magic condition 
for order 0: 



A { °\Y) = kA(Y) (D2) 

with 

4°Vy) = - f A(t,Y)dt. (D3) 
T Jo 

And for order 1, we get 

dti [ 1 dt 2 A(ti,^(t9,y)y)) - A^(A^(Y)Y)+A(t 1 ,Y))A(t 2 ,Y) - A^(Y)A^(Y) = 0. (D4) 
o Jo 

The expansion is closely related to the magic conditions obtained by the Magnus expansion, but extra terms are 
present due to the non linearity of the equation. To stress out this similarity, we will use these non linear magic 
conditions on the mean-field dipolar Hamiltonian. It is indeed well known that the mean-field dipolar Hamiltonian 

t 

with Bdi the local dipolar field which is 

used with quasi classical density matrices p = ®iPi is a formalism directly equivalent to the classical formalisrrP^fil. 
Indeed the magic condition at the order zero in this case is 

= \ [ d*i E bi/TiA'Wfi ) (D7) 



We thus directly see the relation with the quantum condition and find that the geometrical magic condition (27) is 
sufficient to ensure that the classical condition 

Hf ] - k,H d (D8) 

is verified. It is trivial to see that the geometrical magic condition is indeed equivalent to the classic magic condition 
at order 0. 

For order 1, the situation is more complexe as we have 

J ° Jo i7 tj k=ti (D9) 

+ t f i A\t 1 )(f j )%A'(t 2 )(f l ) - ^A'W^y^A'^^) 
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To obtain this expression we use the fact that idt(I) = ([H' d (t), I]}. The expression obtained in the quantum 



case is the same as (D9| without the brakets. We thus see that the geometrical magic criterion of order 1 is a 
sufficient condition to verify the magic conditions in the classical case at the order 1. In appendix [EJ we proove 
that these conditions are in fact equivalent. Magic conditions of higher order can be derived the same way, the only 
differences between quantum and classical conditions beeing those averaging operators appearing in the classic case. 
Those operators having no influence on the geometrical tensor, we deduce that the geometrical magic conditions are 
always sufficient condition both for classic and quantum magic conditions. The equivalence of classical, quantum and 
geometrical conditions at orders greater than one is still an open question. 

Appendix E: Proof of the equivalence between the geometrical magic conditions and the classical magic 

conditions 

To demonstrate the equivalence between the magic geometrical conditions and the classic geometrical conditions at 



the order one, we start from eq. ( D9 ) . Let us first simplify this condition by remarking that this condition is fulfilled 



if it is valid for any vector on the Hilbert space. Hence we can use a braket on condition (D9| to get an equivalent 
condition. Now let us remark the following relationship 



/ dii / 1 di 2 A^A'j") = / dh f 1 dt 2 {A' ab (*i)A'cd(t 2 ) + A' ab (< 2 )A' cd (t 1 )} (El) 
Jo Jo Jo Jo 



so that we can define the tensor T: 



T a ,b,cd = / d*i f 1 di 2 {A / ab (t 1 )A' cd (t 2 ) - A' ab (i 2 )A' cd (i 1 )} (E2) 
Jo Jo 



This expression allows simplifications in (D9| and we get the following expression for the classical magic condition 
at the order one 

o= E E 6 «E 6 « Ta > b . c . d {<4a)([4b ) ^,d])(4c) + (4a4c)(4b)<4d)}- (E3) 

(a,b,c,d)G{x,y,z} i^j k=£l 

We simplify the first part by noting that 

(4,b=4d]) = (4e)£bdefy- (E4) 

In the contest of classical mechanics, the test vectors are limited to separable states which have the property that 

(4a4,c> = <4a)<4,c) (E5) 

if i ^ h. This symmetry implies that due to the structure of T, the term associated with this case is always null. 
If i = k we have 

(Ii,alk,c) — (Ii,e}£ace8ik- (E6) 

The classical conditions can thus be rewritten 

°= E E 

hjbkl Ta, b ,c,d {£bde8jl(h,ii)(Ij,e)(Ik,c) + £a.ce8ik(Ii,e)(Ij,b)(Il,d)} (E7) 

(a, b ,c, d )G{x.y,z} 

Simple permutation of indexes using the symmetry bij = bji allows to write the condition in the form 

0= E h » h 0l E Qa,e^(4a)<4e)(4d) ( E8 ) 
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where Q a ,e,d = Ja^.c.dEbcd which has the property 

Qa,e,d Qd,e,a- 



(E9) 



The equation (E8 1 gives the classical magic conditions of order one while the equation Q a .e.d = gives the geometrical 



magic conditions. It is trivial to see that the geometrical conditions imply that the classic conditions is verified. The 
opposite is more complex but feasable by supposing a sample of 3 spins with the different separable states. 

Let's start with a sample prepared in the state \z + , z + , z + ) (the order n in braket notations corresponds to the spin 
number n). The classic condition in this particular case is 



= ^ hjbjiQz 



(E10) 



If the positions of the spins are correctly chosen the spatial part is non zero which imply that Q z ,z.z = 0. Conse- 



quently by changing the test vector we also get two other geometrical conditions: Q x ,x,x = Q 



y,y,y 



= 0. 



We now use a sample prepared in the state \x+, x+, y+) and compute the classic condition. After some algebra we 



get 



= ai<3 X:y!X + a 2 (Q 



with 



Qy,x,x) + o^Qy x y + conditions already verified 



(Ell) 



a 2 — 

«3 = 



26i 3 6 23 + 6? 3 + 6l 3 
612623 + 612613 



(E12) 



H 3 . 



This condition must be met for any choice of position. The difference in the change of ot\ and a 2 with the choice of 
the spin positions imply that the previous condition is valid for any choice of a\ and a 2 ■ We deduce the geometrical 
conditions Q x ,y.x = and Q x ,x.y = (remembering that the sufficient geometrical conditions has the symmetry 
(E9)). By choosing other directions in the test vector we get the geometrical conditions for every choice in the tensor 
Q with 2 directions. 

Finally, we choose the test vector |x+,y+, z + ) and get the condition: 



= /?l(Qx,y,z + 



y,z,x 



'x,z,y 



) + &(Q 



z,x,y 



with 



Qy,x,z) + conditions already verified 



(E13) 



A= 612623 

P2= 613623 

h=bi 2 b 13 . 



(E14) 



The same reasoning applies here: The position of the 3 spins can be chosen arbitrarily so that the coefficients f3 are 
independent. We know already that the conditions with less than 3 coordinates must be canceled in order to have 
a magic sequence so that we get finally all geometrical conditions (using the symmetry of equation (E9 1: Q a .e,c = 
for any choice of a, e and c. We thus proved the equivalence between the geometrical and classical conditions. The 
equivalence with the quantum conditions is likewise proovable. Extensive analysis using Mathematica show that they 
are only 6 independent first order geometrical magic conditions. By using eq. ( 24 ), it is easy to show that quantum and 



classical conditions can be written as linear set of equations involving the geometrical conditions. LU decomposition 
allows to show that the number of independent geometrical conditions is 6 both in quantum and classical case. In 
the general case of order n > 2 (a case not considered in practice in this paper), the equivalence between classical, 
quantum and geometrical magic conditions still remain to be proven. 



* Electronic address: emmanuel.baudin@lpn.cnrs.fr 
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We recall that we considered h = 1. 

We take the convention that Oi is the rf field amplitude before the rotating wave approximation. 
See appendix [b] for the occupation density explanation. 

In all the presented simulations, the maps of the magnetization distribution confirm that magnetization is indeed mainly 
homogeneous around the instant of the echo. 

Comparisons of magnetization maps between ideal and disturbed evolutions confirm quantitatively this point. 
The equivalence presented below can be found using the Bloch-redfield equation directly with the Picard iteration scheme. 
The mean-field formalism is presented here in order to evidence clearly the similarities between the quantum and classical 
conditions. A presentation which is easier as quantum and mean-field notations coincide, whereas classical notations are 
different. 



